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0 15 ۷۵۲ 015 0۶ 6و0مام 1 .۱65۵2۲0005 ترصقوط رها فصمل موه فقط معبمآتم‌تبهع طا داعم عصقام هم ممتاهتاله۴۱۲۷ 
کنط 1 .۲6۵106 همع ما دومن ۵۶ ۲۵5۵0886 ظ) ع‌صتامن0ه تم 4ص عمتصصعان م۲ مود مها میهع200۳0۵ عصا مصتصماعل 
۲ ۵۱20 ۵۴ مه هه ممصمصصمگه مفطا رم فتمامه قمع ۵۴ ]مع]له مطا «ته م) ماماتوومج 11 فص ود ما2 
0 ,۷7۳۵۵6 15 صو1۳ ما امتالمين آمتنطانمز نموه اصم) تم مرا دمص مط ممصلو. .رمصموط قصع مصصت و15 عصل‌مممی رها ممتام‌حصاقومم 
معط ما معتطلنه‌ننيه ماطفمتهاعته مه وتتهعی 1000 صا عام۲ اصم)تممرصصا صج مق فامی همه ۵ اممجصم‌مهمفظ تمم۳۲۵ 
عط مع2عععصا ما حمتاعج رصح 0ص روامنال0۳0 ما همه )۵2 و1 مق رز ماممعج ما 10۲ 1000 ۵۲ ممتا5۵ صتقصه م1 .تمتاصنامم 
معط[ رم ره »رعآم‌صمی 2 ود ممتاهژمفصه ۳۲۷2۵۵6۲ .وعمتتامکع۲ 5011 24 ۳۷26۵۲ 60تحصناً ما فمل وتحووع۵6ظ و عمط ۵۶ 71۵10 
0ص ۵و۲ رحمتامن124 بل‌عموو 4ص ,تاتکتصنط ,تمرم وق طمتای وعماعع؟ متامصصتای فتامته7 م۵ ول‌صممرعل 24 2۳۵6895 
,0227۷10) بفقک رقتصنا 0۶ ماقل ۵۱60۲0108162 انم عصتفن ود رول‌بتاه کصعومزم مط ا ۵۲6۵10۲۵ .همع کصعاح ؟ه معهاه 
-۳۵۵۲۵- ۳۸۵۵ 6 0۶ وااناعع۲ مط) وم 260 عماله۷ متام تلم قمع امموجنن «انعل معهع7ه مطا رفعمتاهای ۷220 مه و۷۲2 
0ص تصهصه-و۲627۵ع ۲۱2۲) 0مطامجظ متنه تمه و0مطامصه وه معط ۵ برمج‌تنتاهعع فص 0ج 2006160 مه 0مطامجط طتممم]۷( 
۳۵۵-6 ۷۴/16۲ 60۳۲۵2۲60 ۳/۵۲۵ (عاص۱2۱ صرح ع۲ . ,ماو ۲ -0اعمز:ظ۳) وهمطاعصظ ممت۲2012 مععطا مج رماللژین-مصهاظ 
1۰ 10۲ 


۵5 2:0 عاه1۷]۵)6۳1 


ممناهتصلی صز امد ومتس)هتامظ مط) ۵۶ ماه مه رمحتناهمع فا مققتالجبه م6 0م)6۵8000 ۷۵۵ وتو جمومتج مط 1 
56 1 20100امرکصها0 ده مممم۲۵۲۵۲ ممتاقصتایی تم فل‌مطامصه کصمتعن مصتویا روعقعومتها همه ممتاهتامکصهتاممهن 0۴ 
-تصصعو صرح 0تحصباط باه تصصهد بتنه) فمافصصتاه کصمعتل تاو (0هصطمه۱۷ جع ۷۲۵20 ر)حاعق! بص۵221 بقتصتا) مصمتاهاه 
8 3 متقتاده ما 880 فممتاقتوه م1 .ممتامتلمس عمط رم صق. صا ۵0ع0 همع 7۵۲۵ (0تصصناط 
0 ,(۳.۲) 2110۲ -۲۱۱۵۹۲۱۵۷ ,(0.ظ) 011001)-۱2207ظ و(.۲۱) تصقصهو-وع7 ۵۲۵2 ۲۱2۲ :مه ولباو متط) ما ممتاه تام کصهام مهم 
۲ ۷/۵۵۵۶ 10۲ وومتاهاه 0مصماصمجطر فص آم هه فطع هد 0عتهم‌صم م۵۲ وااناوع1 فطع مصفط ,(۷۲۸16) ۱۷۲21 40صه (1) 
-طامه۱۲ 0صج (۱۳۵/۲6۴) عم عتعبوو معط )۲00 تمه بر کصمته‌تگلهمی ممتعصماه :ع مت امعم متمعماتن متا فتاهاه تمه 
۰( 10067 6116 


1650166 200 


وز ماع 0016زین- معا عمط مر ممقصمئط 4صه ۴۲ ممتدامصته ۵۶ رخ )صعئه‌تللهمی صمتاعمداوه عط که عسلعه مط] 
ام ۲۵۲ وماه۷ منام)ن معا 40صه ۲۳۴۵۷۲9۳۲ م1 .026 مطا موه ممتاقآمتنمی ۵0مع 2 و7«مطه طمتط۲۳ بعقصی مه م61 
مص دماهتصا منم ۲۵5۵6۵ موم ها عومآم صع 0-20 ۵۲ مومع مطا صا وی فصمتعاه 1176 معط فصح وتمتامصصهتهم 
-۳۱2۲۵۲6۵۷6۵ عط مز ۲2 ۶ه مبلو* عطا ر‌صقط تعطان قطا هم .ععقصمنطا فصه ۲۲ متعایسصته ما وتانطه و اعمص ومتمناوخ۸ر 
,0.5 صقطا مرمحط و مها میگنآم ناو -طوعا فص 0-10 ۲ معصه مطا مر ۱۱۴۳۱9۳ یعقوم م6 عومام وعقصمتها ۲۵۶ 0مطامحط تصعه 
0 و1 ۷۷۳۵۵۵۲ وفهصومزها هه ۲۲ ۵۶ ممتاهتله7ه معط ی 1۵06 ۱۱۳۴۵۷۲۵ 1 .صمتهاناصته 0ممع 2 ومهع01صا منظ۳ 

و1 5و2صصمزها فص ۶0۲ صرح تمه ود ٩1‏ 0۲ تصهصهو-وع7 ۲۱2۲۵۲۵2 ]جمجاج مه مامح لین مصقاظ فطا 
2 20 هه ما لقاع ۵۲ ما مدمه )صمم‌تگتممه ممتاقصهانه ,0-30 ۵۲ معصه: معط صزً ۱۱۳۴۵/۵۳۲ ۳0 00طاعه 1126 1 
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۲ ,۸150 .50۵11005 1۷6 211 ۵1 عععجصمتها 0صنه ۲۲ میقاناصصته م) ۵0وتا 9 صقع عصم هن عوم01 ۵۲ مهم ۵۶ 1۳06 عکدآماناو-طوها 
۰ 11۷۶ 211 طا فمنیله۷ لمماو20]٩‏ عاصامم0عع2 2۷۵ظ ع0مطامصظ کلن۱۸1 صرح تمانوه 1 م1 ومز۳ظ رصمتامانصته ععقصمته 
,005طاعحط حمتاهن 1۵ 10۲ )عمط وعمحصمتوا فصج ۱ ممتامصتادی ۵۶ (ل) ]صمئم‌تلهم ممتاحصدامونه که ملع مطا جه ۳2860 
ات۱۷۵ ۶ رم0صا ۱۱۳۴9 ۵۲ معمامم۳۵۲ .طونط و1 و0مطامصه ممتاهت0ع مععط الم صز 0242 مط صمعاعط ممتلهاه من مطا 
ص صرح 2000 ود عطقم مه هتصرا صا فص معمتهه ممم ود ممتاهای ۵221۵ ما ممتاهتاله7ه ۲1 ص1 قمع ۶۵ 0مصاع2 
(۷۲۵70 200 ۱۲۳۴۵12 ,۱۷28820 رصاعه رصذ0۵227) عصمتعاه 11۷6 211 ها عفحصصمتها اتامجاق 20 2006۲266 و ۷22۵0 0ص2 1۷28۳20 
جمتاهای ۷22۵0 1 ممتاهت اه ۲۱۱ صا عمط 10 0مطامحظ ممارره ]1 -ماعم۲۳۱ ۵۲ ممدامهمتهم تمه مطا رروامععاهه اصماممره 15 
بلمططعه]۱۷ رکطعف بص022) معمتامء م7 اه ما کصمااهعه وا ععقصهمنما غتامحا2 رتممه ود فصمتاهاه عطا ۵۶ اوع۲ مط 0هصج 2۵00 15 
عصمتاهاه ۱۷۲2۵20 0ص )امه رهتصا ما متاهتا هه ]۲ صا ۳/۵۵۵۲ ۲0۲ 00دامجور م0۲ 1 0۲ ۲2۵۸۵ 6۲0۲ ۰1۳8۵ ۷۵220(۰ 2۳00 ملصصرت] 
بلمططه۱۷]2 روف رصن ۵22)) معومتامء 1*76 21 ما کصمااممن و1 وققصصمتها مامماه صرح 00ج و 220 ۷ 20 022271 ما صرح 2000 15 
۷۵70(۰ 200 ۲12 
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عط صز ۱۱۳۵/۵۲ رعصه ما عفملم مسله ۲ ده ماع عا0نین-مصفاظ همه 0عصنمهه اوعد معط ما تلع 
و010 ملح۷» ۲۲ از 00حطاع م۲ مد عمجم ما عومام عها ورتلمایی- مها 0صه ر8ممع ما اصع1[6عع۵) 0-2066 0۲ ۲226 
7 2000 2 5۳07/۵0 ۷۷25 088 10 01056 »00ص من وود 0ص اصعمااهعه۵) 0-1056 0۲ ۲226 مط) م1 ۱۱۳۴۵۵۲ ,عم 
۶ حمتاحصتایی ۵۲ و0مطامصه معمطا از ففهصهمتها مه ممتاهم‌کمه‌تامممه ۵ ممتاهامصله ما افص ومتت)قتامخظ 0۶ 
۰ مهن ۸۵ ۲60همررمی منامام کصهتام 72و 
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چکیده 


ارزیابی مدل‌های گیاهی در بخش ؟ شاورزی تو سط ب سیاری از پژوه شگران انجام شده | ست. تعیین مدل گیاهی منا سب برای برنامه‌ریزی و 
پیش‌بینی واکنش گیاهان زراعی در مناطق مختلف ضروری است. این عمل سبب می شود با صرف هزینه و وقت کمتر بتوان اثر عوامل مختلف را بر 
عملکرد و کارایی مصرف آب گیاهان بررسی کرد. با توجه به اينکه معادله ۸۵۸6-56 بعنوان روش مرجع برای برآورد 1 در مدل ۸00۵6۲00 استفاده 
می شود و به دلیل تعداد ورودی زیاد | ستفاده از آن د شوار | ست. روش‌های دیگری همچون روش‌های دمایی و 5 شع شعی وجود دارد که با حداقل داده 
ورودی می‌توان ]17 را با همان میزان دقت برآورد کرد. با توجه به اهمیت این موضوع» تحقیق حاضر به منظور بررسی دقت و کارایی مدل ۸۵60۳ 
در شبیه‌سازی تبخیروتعرق و زیست‌وده. تحت تأثیر روش‌های مختلف دمایی (بلانی-کریدل و هارگریوز-سامانی) و تشسعشعی (پریستلی-تیلور 
مککینک و تورک) برآورد تبخیروتعرق مرجع در پنج اي ستگاه (ارومیه» قزوین» ر شت. یزد و م شهد) و چهار اقلیم (خ شک نیمه خشک مرطوب و نیمه 
مرطوب) مختلف در ایران و برای گیاه گندم انجام شد. طبق نتایج» روش بلانی-کریدل با مقدار 157 بیشتر از ۰/۵» 1150۷1512 در محدوده ۰-۱۰ در صد 
(عالی) و شاخص 8 نزدیک به یک (۰/۹۹) و روش تورک با مقدار 1*7 بیشتر از ۰۰/۵ ۷12۷18۳ در محدوده ۱۰-۵۰ درصد و شاخص 8 برابر با 
٩‏ روش‌های منا سب برای شبیه سازی تبخیروتعرق در تمام ایستگاه‌ها بودند. در مورد شبیه سازی زی ست‌توده» روش‌های بلانی -کریدل و هارگریوز- 
سامانی با مقادیر 187 برابر با ۰۰/4 ۱11152 در محدوده ۰-۱۰ درصد (عالی) و شاخص ۱8[ برابر با ۰/۹۹ بعنوان روش‌های دمایی مناسب و روش‌های 
پریستلی-تیلو, مک‌کینک و تورک با آماره‌های **1 برابر ۰/۹ ۱۷18۷5 در محدوده ۰-۱۰ درصد (عالی) و شاخص 78 برابر با ۰/۹٩‏ بعنوان روش 
های شءشعی منا سب انتخاب شدند. در این پژوهش دقت خوب مدل ۸0۵600 در شبیه سازی تبخیروتعرق و زیست‌توده با این روش‌های برآورد 
تبخیروتعرق نسبت به سایر روش‌ها نشان داده شد. 


واژه‌های کلیدی: بلانی-کریدل, تبخیروتعرق مرجع. تورک. مدل گیاهی, ۸۵6:00 
مقدمه زیرزمینی می شود. ضمن آن که تخمین کمتر از حد موردنیاز گیاه نیز 
باعث اعمال تنش رطوبتی کنترل‌شده به گیاه شده و در نتیجه کاهش 
مح صول را به همراه خواهد دا شت. با توجه به اهمیت تبخیروتعرق و 
کاربردهای متنوع آن در علوم مختلف محاسبه مقدار آن به ویژه 
تبخیروتعرق پتانسیل و واقعی از اهمیت بالایی برخوردار است. 
تبخیروتعرق در کل جهان حدود ۶۰ درصد از بارندگی زمین را به 


تبخیروتعرق از عمده‌ترین اجزاء چرخه هیدرولوژیکی است که 
تخمین درست آن در طراحی و مدیریت سیستم‌های آبیاری» مطالعات 
منابع آب و موارد مشابه دیگر از اهمیت زیادی برخوردار است. تخمین 
بیش از حد آب موردنیاز گیاه ضمن هدر دادن آب آبیاری باعث ماندابی 
شدن ارا ضی» ش ستشوی موادغذایی خاک و آلوده نمودن منابع آب 
۱ و ۲- به‌ترتیب دانشیار و دانشجوی دکتری» گروه علوم و مهندسی آب دانشکده کشاورزی و منابع طبیعی» دانشگاه بین‌المللی امام خمینی(ره» قزوین» ایران 


(*- نویسنده مسئول: ۴ .موه ه) 1۵62۵1 :اتمظظ) 
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جو برمی‌گرداند (2006 ,1026 200 01). بطور کلی معادلات 
محاسبه تبخیروتعرق براساس تابش خالص خورشید (,48 حداقل و 
حداکثر دماء سرعت باد و رطوبت نسبی هستند ( ,.۵1 21 ۸۵116 
8 هرگونه تغییر در پاامترهای آقلیمی ورودی موجب تغبیر در 
می‌شوند (2014 .۵ 61 1۱6۵016-561۲280 ۰۷ ۷۷۷۲ و همکاران 
(2019 .۸1 ۶۶ ۷۷۷) تغییرات مکانی و زمانی در نیاز آب آبیاری گندم 
را با اسستفاده از بارش موثر و تبخیروتعرق محصول مورد تجزیه و 
تحلیل قرار دادند. آنها نیاز آب آبیاری را محاسبه کرده و دریافتند که با 
نیاز آبی محصول متفاوت است و به شدت به در د سترس بودن داده 
های هواشناسی بستگی دارد. در زمینه محاسبه تبخیروتعرق و تاثیری 
که روش‌های مختلف محاسبه آن بر عملکرده زیست‌توده و نیاز خالص 
آبیاری مح صول دار مطالعات زیادی در سطح جهان و ایران انجام 
گرفته است. براساس نتایج پژوهشی که صفری و همکاران ( 521271 
1 ,. 61) به منظور وا سنجی و | صلاح ضرایب معادلات برآورد 
تبخیروتعرق بر روی چهار ایستگاه سینوپتیک ایران با اقلیم‌های 
خشک» نیمه خشک» مرطوب و نیمه‌مرطوب انجام دادنده روش بلانی- 
کریدل بعنوان روش برتر برای محا سبه ]17 در هر چهار اقلیم انتخاب 
شد. در پژوهشی دیگر, سه روش (بلانی -کریدل» تورنت وایت و 
هارگریوز-سامانی) برای برآورد 171 پتانسیل در شهرستان امیدیه مورد 
مقایسه قرار گرفتند که بر اساس نتایج روش بلانی-کریدل ذسسیت ید 
سایر روش‌ها با بی شترین ضریب هم ستگی و کمترین میزان خطای 
استاندارد با روش ۳۸۵0-56 از دقت بیشتری برای برآورد تبخیروتعرق 
پتانسیل برخوردار بود (2014 ,1۵۷00۵ 20 طعتععش). نتایج حاصل 
از مطالعه‌ای که در ایستگاه سینویتیک زاهدان انحام شد نشان داد 
روش‌های بر پایه نتقال جرم از نظر آماری ضعیف‌ترین عملکرد را 
ذسبت به سایر روش‌ها در مقایسه با روش استاندارد ۳۸0-56 نشان 
دادند؛ اما مدل‌های دمایی و تشء‌شعی مانند مدل تورک» جذسن-هیز, 
هارگریوز و بلانی-کر یدل می‌توانند جایگزین مناسبی برای روش 
ترکییی نسیتاپیچیده ۳۵0-56 برای اقلیم گرم و شک شهرستان 
زاهدان باشند (2017 ,ح02025ع0120ه>). در پژوهشی که در 
ترکیه انجام شد. روش هارگریوز-سامانی برای این منطقه اصلاح شد. 
آنها میزان تبخیروتعرق را برای ۲۷۵ ایستگاه هواشناسی و هفت 
منطقه جغرافیایی با استفاده از روش های ۳۸۵-56 و هارگریوز- 
سامانی به د ست آوردنده سپس ضرایب معادله هارگریوز- سامانی را 
واسنجی و معادله واسنجی شده سرعت باد را به‌عنوان یک پارامتر 
ورودی اضافه کردند. نتایج پژوهش آنها نشان داد که روش هارگریوز- 
سامانی اصلاح شده تخمین خوبی نسبت به داده‌های واقعی لایسیمتر 
دارد (2017 .اه 61 تعصدطاهن) 

مهمترین محصول کشاورزی در ایران گندم است و به همین 
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دلیل مدیریت صحیح مزارع گندم. نقش مهمی در امنیت غذایی و 
کشاورزی پایدار در کشور دارد. منبع اصلی غذای مردم در ایران» گندم 
و فرآورده‌های آن است و هر اقدامی برای بالابردن عملکرد محصول 
گندم با توجه به منابع محدود آب و خاک ضروری است ( 6 ۸۵0201 
7 ,.0). مدل‌های گیاهی در جهت بهینه سازی عوامل مدیریتی» 
اجازه شبیه سازی مح صول را می‌دهند و قابلیت کارکرد سناریوهای 
مختلف را دارند. یکی از مدل‌های گیاهی پر کاربرده مدل آکواکراپ 
است که برای محدوده وسیعی از محصولات زراعی شامل محصولات 
علوفه‌ای. سبزیجات» غلات میوه‌هاء م<صولات روغنی و غده‌ای می 
تواند به‌ کار گرفته شود (2014 ,0/06 ۷). مدل آکو کراپ در 
تحقیقات مختلفی در مناطق مختلف در شبیه سازی محصول گندم به 
کار رفته است. 

اقبال و همکاران (2014 ,.41 ۶ 10021) مدل آکواکراپ نسخه 
۱ را برای محصول زمستانه گندم در دشت شمالی چین واسنجی و 
صحت سنجی کردند. نتایج بررسی آنها نشان داد عملکرد زیست‌توده 
در شرایط مختلف کم آبیاری با دقت منا سبی تو سط مدل برآورد می 
شود. همچنین این مدل تو سط مابهودی و همکاران ( نط۷]2020 
4 ,.21 اع) در منطقه نیمه خشک آفریقای جنوبی برای یکی از 
گیاهان گرم‌سیری آن بنام تارو مورد ارزیابی قرار گرفت و زیست‌توده 
با ضریب تعیین ۰/۹۹ و ريشه میانگین مربعات خطا برابر ۱/۷۴ تن در 
هکتار برآورد شد. مدل آکواکراپ برای تعدادی محصول و برخی از 
مناطق ایران نیز موردارزیابی قرار گرفته اسست. این مدل در منطقه 
کرج برای پیش‌بینی عملکرد گندم و سویا در شرایط کم آبیاری نتایج 
قابل قبولی ارائه داد ( 24 2022۵06ظ :2010 ,.۵ ۶۶ طع220خر 
2 ,1۵0۳121 2721). 

زندیار سا و همکاران (2016 ,.۵1 61 22001-۳2752)» در پژوهشی 
مقادیر رطوبت خاک در عمق‌های مختلف» ماده خشک و عملکرد دانه 
گندم زمستانه در پنج تیمار آبیاری دیم» ۰۰/۵ ۱۰۰/۸ و ۱/۲ برابر 
آبیاری کامل را طی دو سال زراعی ۱۳۸۸-۸۹ و ۱۳۸۹-۹۰ واقع در 
اراضی دانشکده کشاورزی دانشگاه شیراز توسط مدل‌های آکوا کراپ و 
" شبیه سازی کردند. داده‌های اندازه‌گیری شده در سال اول 
زراعی برای وا سنجی مدل‌ها و در سال دوم برای اعتبار سنجی مورد 
استفاده قرار گرفتند. صحت برآورد رطوبت خاک توسط مدل‌هاء مبنای 
مقایسه دقت برآورد مقادیر شبیه سازی شده تبخیروتعرق واقعی گیاه 
قرار گرفت. بر | ساس نتایج» در مدل ۷۷۹۱۷۲ مقدار رطوبت خاک در 
لایه‌های مختلف عمق ريشه در سال اعتبارسنجی با مقدار نرمال شده 
ريشه میانگین مربع خطا (۱18۷1517) برابر ۰/۱۴ خوب برآورد شسده 
است ولی مدل آکواکراپ در برآورد آن د قت کمتری داشت 
(۱۱۳۷5۳-0.26). مقادیر برآورد تبخیروتعرق گیاهی در تیمار آبیاری 


رمضانی اعتدالی و صفری. رمضانی اعتدالی و صفری. ارزیابی تاثیر روش‌های مختلف برآورد 1710 در شبیه‌سازی تبخیروتعرق واقعی .... ۵ 


کامل» در هر دو مدل نزدیک به هم ولی با فزایش تنش رطویتی» از 
دقت برآورد مدل آکواکراپ کاسته شسد. مدل ۷۷[۷ مقدار ماده 
خشک و عملکرد دانه را با 21[۷]8۳ به ترتیب برابر 9۰/۱۵ ۰/۱۸ 
خوب برآورد کرده است و مدل آکوا کراپ با ۷1517( به ترتیب برابر 
۹ ۰/۳۹در برآورد آنها دقت کمتری داشت. 

امیری و همکاران (2015 ,.1 6 نتنتش) توانایی مدل آکواکراپ 
را برای پیش‌بینی عملکرد دانه و زیست‌توده کل گندم زم‌ستانه تحت 
تنش آبی؛ به مدت دو سال برر سی کردند. تتایج ذ شان داد مدل در 
پیش‌بینی عملکرد دانه و زیست‌توده محصول دارای دقت قابل قبولی 
است. عبداله‌زاده و همکاران (2019 ,.۵1 ۶۶ ۸00۵1122206۳) میزان 
تبخیروتعرق واقعی و نیاز خالس آبیاری محصولات گندم. جو و ذرت 
دانه‌ای را در د شت مغان با | ستفاده از مدل آکوآکراپ تعیین و نتایج را 
با مدل کراپ وات و نت وات مقایسه کردند. براساس نتایج» 
تبخیروتعرق و نیاز آبی حاصل از مدل آکوآکراپ در مقایسه با کراپ 
وات کمتر و نتایج در مقایسه با نت وات برای گندم و ذرت دانه‌ای 
مقادیر بیشتر و برای جو کمتر برآورد ده است. در پژوهشی که 
جرعه‌نوش و همکاران (2019 ,.۵ ۶ 1076015۳) برای شبیه‌سازی 
عملکرد گندم و تعیین تاریخ کشت در | ستان فارس تو سط آکواکراپ 
انجام دادند» نتایج نشان داد که مدل دقت بالایی در سبیه‌سازی 
پوشش تاجی» زیست‌توده و عملکرد دانه گندم دارد. نتایج مطالعه‌ای 
در دهلی هندو ستان نشان داد که مدل آکواکراپ دقت قابل قبولی در 
شبیه سازی عملکرد دانهه زي ست‌توده و کارایی م صرف آب در ارقام 
مقاوم و غیرمقاوم به شوری گندم دارد. در این تحقیق م شخص شد 
که تونایی مدل در شبیه سازی عملکرد بیش از دو پرامتر دیگر است 
(2014 ,.۵1 6۶ تقصناک) 

پژوهشگران بسیاری در سراسر دنیا روش فائو -پنمن- مانتیث را 
در مقایسه با لایسیمتر به عنوان روشی دقیق معتبر و استاندارد برای 
محاسبه تبخیروتعرق گیاه مرجع معرفی کرده‌اند که در این پژوهش 
نیز به عنوان روش مبنا مورد استفاده قرار گرفته است ( 6 16۱027 
:998 ,۵1۱72060 0ص ۲۵۸067 2002 .601 تبخیروتعرق فرآیندی 


پیچیده و غیرخطی است و به عوامل متعدد اقلیمی نظیر دما رطوبت» 
سرعت باده تابش, نوع و مرحله رشد گیاه و غیره وابسته اسست 
(2009 ,.۵1 6 82(201). بنابراین» در تحقیق حاضر با استفاده از 
داده‌های هواشناسی روزانه مربوط به ایستگاه‌های ارومیه» رشت» 
قزوین» مشهد و یزد مقادیر تبخیروتعرق متوسط روزانه بر اساس 
نتایج حاصل از روش فائو-پنمن-مانتیث مدل‌سازی و دقت دو روش 
دمایی (هارگریوز-سامانی و بلانی-کریدل) و سه روش تشعشعی 
(پریستلی-تیلوره تورک و مک‌کینک) با ۲۸0-56 برای گندم مقایسه 
شد. بنابراین می‌توان نوآوری این تحقیق را ارزیابی اين مدل در برآورد 
زیست توده و تبخیروتعرق واقعی را با استفاده از سایر روش‌های برآورد 
تبخیروتعرق که گاه متفیر های اقلیمی کمتری ذیاز دار ند در مقا بل 
روش پایه ۲۸0-56 که در مدل بعنوان روش اصلی و پایه شناخته 
می‌شود بیان کرد. 


مواد و روش‌ها 
مناطق مورد مطالعه 

ایران با وسعتی حدود ۱۶۳۸۰۰۰ کیلومتر مربع و متوسط بارش 
سالانه ۲۴۱ میلیمتر در نیمکره شمالی» در قاره آسیا و در قسمت غربی 
فلات ایران واقع شده‌است (2004 ,.۵1 61 0211271)). این کشور بین 
دو نصف‌النهار ۴۴ و ۶۴ شرقی و دو مدار ۲۵ و ۴۳۰ شمالی قرار 
گرفته است. حدود ٩۴/۸‏ در صد از سطح آن, در زمره مناطق خشک و 
نیمه‌خشک با ریزش های جوی کم و تبخیروتعرق ز یاد قراردارد 
(2004 ,.]۵ 61 تجدل0)). برای انجام این پژوهش. ایران را براساس 
اقلیم نمای کوپن به ۴ اقلیم مختلف (خشک» مرطوب. نیمه‌خشک و 
نیمه مرطوب) دسته‌بندی و به‌ترتیب شهرهای یزده رشت مشهد» 
قزوین و ارومیه بعنوان نماینده این اقليم‌ها انتخاب (شکل ۱) و آمار 
هواش ناسی آنها مورد استفاده قرار گرفت 
(۰1۳/]2/08605/156320ع05://۵216.۰0۳). خصو صیات ایستگاه‌های 
هواشناسی در جدول ۱ ارائه شده‌اند. 


جدول ۱- مشخصات ایستگاه‌های هواشناسی مورد مطالعه 
مصماهای ۱۵۲۵۵۲۵۱۵۵ 0۲ وانهاع( -1 12016 


ارتفاع ار سطح دریا عرض جغرافیایی طول جغرافیایی 
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جدول ۲- معادلات مورد استفاده در پژوهش 


12016 2- ۲۳/۵۳۵۵۵0۵5 ۱۹60 1 ۵ 


منابع فرم معادله نام معادله شماره ردیف 
۱۹ حصم] ممتادن ۳ .۳ 0۲ عصولا و۱۱۱ 
990 
(0999) 1۸0-56 و ۱ ۳۸0-6 1 
(را0.34 + ۲۷1+ ۸ 9 
(1985) تصحصوی 20 ۲۱۵۲۵۲۵۵۲6۵ (17.8 ۲ مگ نگ < معج 0.00231 < و۲1 تصحطه ۲۱2۲2۲62۷65-۹ 2 
(1950) 00016 290 ما۳ ((8.13 چ موم ۲)۳)0۰461 +6 < و1 26 - رمصد۱ظ 3 
5 ۸ 
(1972) 131۲10۲ 0صد 16۲)ومز۳ 7 ۳ 7 23 ۳ < و17 07 - ۳۱[69)161 4 
1 ۵ 
(1957) ۷۲۵۱0۱0۱ 2 - ۳ لس 1 << 7[ ص۷۲2( 5 
5 (7- ۸ 
7 + 23.89۳ 
(1961) 10۳ ۳ ۲ 2 - 7 ۳ 6 


(15 + وس 1) 


روش پژوهش 

معادلات مورد استفاده برای برآورد تبخیروتعرق مرجع در این 
پژوهش از دو گروه دمایی شامل: هارگریوز- سامانی (1.5) و بلانی- 
کریدل (60.() و گروه تشعشمی شامل؛ پریستلی-تیلور[2.7 تورک 
(7) و مککینک (121) انتخاب شدند. فرم اصلی معادلات در جدول 
۲ ارائه شده است. 

داده‌های تبخیروتعرق مرجع بعنوان ورودی‌های مدل برای اجرای 
مدل آکوآکراپ موردنیاز است. به همین منظور متفیرهای هوا شناسی 
شامل حداکثر و حداقل دمای هو حداکثر و حداقل رطوبت نسبی هواه 
ساعات آفتابی و سرعت باد در ارتفاع دو متری سطح زمین از ایستگاه 
های ارومیه» رشت قزوین» مشهد و یزد از تاریخ ۱۹۸۰/۱/۱ تا 


۱ جمعآوری و تبخیروتعرق مرجع با استفاده از روش‌های 
مذکور در جدول بالا محاسبه شد. 


به دلیل ناکافی بودن داده‌های لایسیمتری از روش فائو-پنمن- 
ماندیث (۳۵0-56) با تو جه به این که دارای د قت بالا در برآورد 


تبخیروتعرق مرجع اسست. از برآوردهای این روش بعنوان داده 
مشاهده‌ای و مرجع برای صحت‌سنجی استفاده شد. 


معرفی مدل آکوآکراپ 

اساس تخمین عملکرد محصول در مدل آکوآکراپ راب طه 
و00۳6( و 16292 است که در شماره ۳۳ نشریه آبیاری و 
زهکشی سازمان خواربار جهانی (فائو) ارائه شده است. با اعمال 


رمضانی اعتدالی و صفری. رمضانی اعتدالی و صفری. ارزیابی تاثیر روش‌های مختلف برآورد 110 در شبیه‌سازی تبخیروتعرق واقعی ... ۷ 


اصلاحاتی در آن از جمله تفکیک تبخیروتعرق واقعی (7) به تبخیر 
از سطح خاک (:) و تعرق (1) و نیز عملکرد به بیوماس (ظ) و 
شاخص برداشت (]۳1) استنتاجم شده است (2012 ,.۵1 61 5۵65): 
1 ۲ 

)1- (ج‎ <  )۵ - ۷ 

که در آن 28 حداکثر عملکرد» 26 عملکرد واقعی. ۰1 حداکثر 
تبخیر و تعرق» ۳1 : تبخیر و تعرق واقعی و یک فاکتور تناسب بین 
افت نسبی عملکرد و کاهش نسبی تبخیر و تعرق است. با 
پیشرفت‌های علمی و تجربی در روابط گیاه -آب از سال ۱۹۷۹ همراه با 
نیاز مبرم به بهبود بهره‌وری آب فائو نشریه ۳۳ خود را اصلاح و تکمیل 
نموده و بصورت مدل آکوآکراپ ارائه داد. برای محاسبه عملکرد از 
زیست‌توده» مدل آکوآ کراپ از رابطه ۲ استفاده می‌کند ( ,.۵1 67 1265 
12 
۳ 8 ار < ۷ 

که در آن نآ : شساخص برداشت مرجع (طی مرحله بلوغ 
مرجع را تنظیم می‌کند. مقدار زیست‌توده تولیدی بصورت روزانه در 
مدل آکوآکراپ از رابطه زیر بدست می‌آید: 

« 1 

92 ۲ ۲( 

که ۱۷۷۳۳ : بهره‌وری آب مصرفی گیاه» ۹۹ تعرق روزانه 9 7/10 
تبخیروتعرق مرجع روزانه می‌باشند (2022 ,.۵1 1 5265). مهمترین 
داده‌های اقلیمی مورد نیاز مدل عبارتند از: حداقل» حداکثر و ميانگین 
دمای روزانه. تبخیروتعرق گیاه مرجع (50) و بارندگی. مدل از داده 
های حداکثر و حداقل دمای روزانه برای محاسبه درجه روز رشد به 
منظور تعدیل عملکرد زیست توده بر اثر خسارات ناشی از سرماء استفاده 
مورد نیاز برای محاسبه و۳1 از اطلاعات آماری سالهای ۰ ۱۹ 
ایستگاه‌های قزوین» رشت. ارومیه» مشپهد و بزد استفاده و م۴ با 


شاخص‌های آماری 

در این تحقیق نتایج حاصل از سناریوها با داده ایستگاه‌های مذکور 
برای گندم توسط معیارهای آماری خطا شامل آماره‌های ضریب تبیین 
(82» ریشه میانگین مربعات خطای نرمال شده (۱01811517) و شاخص 
نش- ساتکلیف (5.» مقایسه شدند. ضریب تبیین یکی از مهم‌ترین 
معیارهای ارزیابی ارتباط میان دو متغیر * و ۱1 ست که به صورت بی 
بعد نمایش داده می‌شود. این ضریب ارتباط مسستقیمی با ضصریب 
همبستگی دارد. به این ترتیب که با جذر گرفتن از ضریب تبیین 
می‌توان ضریب همب ستگی میان دو سری مورد برر سی را به د ست 
آورد. همانند ضریب همبستگی هر چه مقدار ضریب تبیین به یک 
نزدیک‌تر با شدء ارتباط قوی‌تری میان دو متغیر وجود دارد. در صورت 
ضرب ضریب تبیین در عدد ۱۰۰ مقدار به دست آمده بیان گر در صد 
واریانس تغییرات متغیر ‏ ست که با متفیر ۷ تو صیف می شود. طبقه 
بندی ضریب پیر سون در جدول ۳ آورده شده است ( ,.01 6۶ 010۲اگ 
7 رابطه ضریب تبیین به شرح زیر است: 

2ج -: 7 
۳ -م 

که ۳1: مقادیر پیش‌بینی شده 1): مقادیر اندازه‌گیری شده و ): 
میانگین مقادیر اندازه‌گیری شده می باشد ( 20 صد/ومححوظ۳ 
89 ,.۳10000000۳1). برای محاسبه ضصریب تبیین از نرم‌افزار 
اکسل استفاده شد. 

شاخص ۷۲5۳۲ سطح تخمین را ذ شان می‌دهد. تق سیم‌بندی 
۲ توسط چمیسون و همکاران (1991 ,.۵1 61 0ععنصصد) 
در جدول ۴ آورده شده است. 


2 -. 1 
7 نصکعت/ 2 ع ۱۲7۲۳68 
10 0 


جدول ۳- تقسیم‌بندی نتایج شبیه‌سازی ضریب پیرسون 
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جدول 4- تقسیم‌بندی نتایج شبیه‌سازی بر اساس :۱۱1۵۱15۲ 
0 02660 کااناعع۲ ممتاهامصته که صملای1]لمعم0 -4 19016 


10-20 20-0 <30 


ضعیف متوسط خوب 


0-10 ,۱( 
عالی نتیجه تخمین 


۶ ان م۱ ...نامع ممتافصتاع۲] 


۸ نشریه آب و خاک. جلد ۰۳۷ شماره ا. فروردین - اردیبهشت ۱۴۰۲ 


ضریب نش-ساتکلیف )٩(‏ یکی از رایج‌ترین شاخص‌هایی است 

که برای ارزیابی کارایی مدل‌های هیدرولوژیکی استفاده می‌ش‌ود 

(1970 ,ع1ذ1هاا5 290 طدعآ)» این شاخص حالتی | ستاندارداز تابع 
حداقل مربعات خطاهاست: 

2«(:ط0) 81 - (صنی) 1۲ ار 


( 05((2) 1-(009) 7ص ابر - 21 ۷5 


دامنه تغییرات این شاخص از م- تا ۱+ بوده و مقدار بهینه این 
شاخص یک است. براساس پژوهش‌های مختلف صورت گرفته در این 
زمینه از جمله مطالعات گاسمن و همکاران ( ,له اه صفصععدن 
7 اگر مقدار ضریب نش-ساتکلیف بالاتر از ۰/۵ باشد مدل 


شبیه‌سازی خوبی داشته است. 


جداول ۵ و ۶ نتایج مربوط به ارزیابی آماری 177و زیست‌توده را 
برای هر ۵ ایستگاه هوا شناسی در سال‌های ۱۹۸۰-۰۰ برای گندم 

براساس جدول ۵ مقدار ضریب تبیین (157) برای شبیه سازی 171 
و زیست‌توده در روش بلانی-کر یدل نزد یک به یک بوده که 
5 نیز در مورد هر دو پارامتر و پنج ایستگاه به ترتیب در محدوده 
۰-۰ و نزدیک به یک است که نان دهنده توانا بودن مدل 
آکوآکراپ برای شبیه سازی ۳ و ژیست‌توده می‌با شد. از طرفی مقدار 
2 در روش هارگریوزسس-امانی برای زیست‌نوده نزدیک به یک 
۱۱۹/۹۵ در محدوده م۱۰ 9 شساخص دا بیش از ۰/۵ بوده که 
ذشان‌دهنده شبیه سازی خوب است (جدول ۵ پا شاخواه و همکاران 
(2014 ,.۵ ۶ ط۳2۵۵12) به بررسی تبخیروتعرق مرجع با روش 
های بلانی-کر یدل» هارگریوز و ترنت وایت برای اقلیم‌های مختلف 
6 پرداختند. نتایج آن‌ها نان داد روش بلانی-کریدل در اقلیم‌های 
مطالعه بیان کرد در شرایطی که در اقلیم‌های مورد مطالعه به دلیل عدم 
دسترسی به داده‌های مورد نیا برآورد تبخیروتعرق مرجع از روش 
۳۸6۵-۵6 ممکن نباشد با بهره‌گیری از معادلات واسنحی‌شده می‌توان 
به برآوردهای مشابه از نظر آماری با روش ۳۸۵6-56 دست یافت. 
همچنین روشن و همکاران (2012 ,.۵ 6 05۳020) در مطالعات 
خود در زمینه‌ی ارزیابی مدل مناسب تبخیروتعرق پتانسیل در ایران» از 
چهار روش تورنت وایت» جنسن- هی بلانی-کریدل و هارگریوز- 


سامانی استفاده کردند و به این نتیحه رسیدند که روش بلانی-کریدل 
همخوانی بهتری با شرایط ایران دارد. 

طبق جدول نتایج گندم (جدول 4۶ روش تورک با ۱151 در 
محدوده سم ضصریب تبیین نزدیک پا برابر با یک 9 شاخص ولا 
یک يا نزدیک به یک برای شبیه سازی ]7 و زي ست‌توده در هر پنج 
ایستگاه قابل استفاده است. همچنین برای شبیه سازی زیست‌توده هر 
سه روش در هر پنج ایستگاه دارای مقادیر آماری قابل قبول است 
(جدول ۶). طی پژوهشی که توکلی و همکاران ( .۵1 6 12721011 
3 در استان خراسان شمالی با اقلیم خشک و نیمه خشک انجام 
دادنده میزان تبخیروتعرق را با روش‌های هارگریوز-س‌امانی» تورک» 
پریستلی-تیلوره مک کینک و تشت تبخیر کلاس ۸ در شرایط کمبود 
های مطرح شده (به دلیل عدم وجود داده‌های ۷+ سیمتری) محا سبه 
کردند. براساس ارزیابی آنهاء این معادلات دقت قابل قبولی در مقایسه 
با مطالعات انجام شده دارند که با نتایج پژوهش حاضر برای اقلیم‌های 
خشک و نیمه‌خشک تطابق (قزوین» یزد و مشهد) دارد. 

محققان با | ستفاده از ۷ روش مقادیر تبخیر-تعرق پتان سیل را در 
نیمه‌مرطوب) هند تعیین کرده و نتایج حاصل را با روش ۳۸۵0-56 
مقایسه کردند. آنها نتیجه گرفتند که در میان روش‌های مورد مطالعه 
روش هارگریوز در دوره‌های زمانی روزانه و ماهانه نتایج نزدیک به 
مرطوب روش تورک (روش تشعشعی) مناسب‌تر است. نتایج حاصل از 
تحلیل عاملی ذشان می‌دهد که سرعت باد مهم‌ترین عامل در مناطق 
خشک و9 ساعات آف تابی مهم ترین عامل در مناطق مرطوب و 
نی مه مرطوب است (2005 ,1600۷۵۲ 20 تنع02صهل). در این 
پژوهش نیز برای اقلیم مرطوب (ایستگاه رشت) روش تورک اولویت 
دارد. 

شکل ۲ و ۳ مقایسه مقادیر زیست‌توده و تبخیروتعرق واقعی 
شبیه سازی شده گندم در روش‌های مختلف برآورد 710 با استفاده از 
مدل آکوآکراپ در ایستگاه‌های مورد مطالعه را نشان می‌دهند. با توجه 
به نمودارها اختلاف زیاد بین مقادیر شبیه سازی در روش‌های مختلف 
با روش پایه ۲۸۵-۵6 وجود ندارد. همچنین براساس آماره ضصریب 
تبیین» همبستگی بین داده ها در تمامی روش ها اعم از د مایی و 
تشعشعی برای گندم وجود دارد. 
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شکل ۲- مقایسه مقدار تبخیروتعرق واقعی شبیه‌سازی شده گندم با روش‌های مختلف برآورد 1710 
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ارزیابی مدل ۸02000 برای گیاهان رایج در یک منطقه نقش 
مهمی در مقایسه عملکرد محصول در شرایط مختلف دارد. گندم یکی 
از محصولات استراتژیک در کشور است. از طرفی, با توجه به آهمیت 
دقت در برآورد تبخیروتعرق مرجع و استفاده از آن برای شبیه‌سازی 
میزان عملکرد و زیست‌توده با نرم‌افزار ۸00۵000 دقت داده‌ها در 
برابر داده‌های مشاهده‌ای بد ست آمده برای هر پنج ایستگاه ارزیابی و 
طبق نتایج حاصله» از بین دو روش دمایی (بلانی-کریدل و هارگریوز- 


محدوده (۲۰-۰) و "12 و ٩‏ نزدیک به یک برای شسبیه‌سازی 
تبخیروتعرق واقعی و استفاده از هر دو روش بلانی-کریدل و 
هارگریوز-سامانی برای زیست‌توده قابل قبول است. همچنین در مورد 
روش‌های تشعشعی» روش تورک برای شبیه‌سازی تبخیروتعرق واقعی 
و تمامی روش‌های تشعشعی مورد بررسی برای شبیه‌سازی زیست‌توده 
مطلوب اند. پس با توجه به پهناور بودن کشور ایران و اقلیم های 
متفاوت آن, پیشنهاد می شود تبخیروتعرق در همه اقلیم‌ها و شهرها با 
روش‌های ارائه شده و يا سایر روش‌های اندازه‌گیری تبخیروتعرق 
محا سبه شود و روش‌های برتر در هر اقلیم م شخص و برای کارهای 


سامانی) برای گندم روش بلانی-کر یدل با میزان ۱1۷/51 در تحقیقاتی معرفی گردند. 
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007 ۱۷۵۲۵6۱ 20 ص12 ,مضه ,م۸۵ م۱۵ ۲۵/۰۵۲/۱۵۵ 2012(۰) ۳,۰ روم‌ز۳۲6 ع ریت1 ,10و۲۱ ,.ظ م,مانالعاو ر.ما رو26 
۰ ۲0۳۵ 

00 ۱۷۵۲۵6۲۱ 24 ص12 ,مظر ,0و م۱۵ ۴۵/۰۵۲/۱۵۵ 2022(۰) ۳,۰ روم۳۲6۲ ع ریت1 ,۲۱8120 ,.ظ ,مانالمای رین رو۴26 
۰ ۲0۳۵6 

۶۲ ۵061۱ 016ا2)زناه ۵ م عمط وی 0صه مصترتلمصط رعصتووع‌عع۸ 2012(۰) ۱۷۲۰ رعتامصه هک[ عک ر.۲ رطعها اهوم دک ریت6 رصحطعومع 
.49-68 :(43)078 رایع «امهو0ع آمعتعر۳۳ ص12 صا ممتاهتاموصهتممديه لهتاصماوم اه ممتاقطدتادم 

(0عاعطاه طعتاع ص ۳/10 صهتوع۳ 

۳۱۵۲۵۲۵۵۷6۵۰ ۵۶ واصعتم11ه0 686 ۱۷۲۵011160 .(2021) ۲۲۰ ,۲۲۵0211 2۵1عصصیم ی رش رتفاهدان صه۸2121 وب رتطه1 162۷ ,.۲ رتتفلهه 
معط ۵۶ ممتاهتم‌وصهتممد۵7 معط عقاو مه فصمتاهنامم اه ۲ مج عاصن۱۷]2 رتمانوه ] -وماممز۲۳ رعا00زین-مصقاظ رتصعصه 
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